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1 Abxtroct 

Three-dimensional supersonic flows of an ideal fluid in the neighborhood of 
bodies formed by being cut out along the streamlines of an axisymmetric 
flow are investigated. The flow consists of a region of isoentropic 
compression and a region of vortex flow. An exact solution with variable 
entropy is used to describe the flow in the vortex region. In the 
continuous flow region an approximate solution is constructed by expanding 
the solution in a series in a small parameter. The effect of the shape of . 
the excission- and the yorticity of the flow on compression of the jet and 
the total pressure loss coefficient is studied. 
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DECELERATION OF A SUPERSONIC FLOW BEHIND A 
CURVED SHOCK WAVE WITH ISOENTROPIC PRECOMPRESSION 

V.G. Dulov , V.A. Shchepanovskiy 
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compression and regions of vort 
variable entropy [1] is used to 
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ee-dimensional supersonic / 72 ♦ 

neighborhood of bodies formed by 
es of an axisymmetric flow. The 
ions: regions of isoentropic 
ex flow. An exact solution with 
describe the flow in the vortex 
region an approximate solution is 
lution in a series according to a 
excision shape and flow vorticity 
1 pressure loss coefficient is 


1. Reference [1] investigates solving axisymmetric equations 
for an ideal gas with variable entropy. This flow is described 
with a common differential equation for characteristic function 
z(h): 


vW r 


’ . [2a 1 • 9 • • 9 / - ' r- ?'• 

2 . Z * 

•«“ Ml -h)~ (a - 1) -y ( 1 + [a.-; 

-.‘A - • -aM* r i)-4“ s ]^'^ 

b = h- (/t (xa* — 4xa + 3a s ) — 4az (a — 1)], 


( 1 . 1 ) 


and total pressure distribution a 
as p o = C(^o - M?)** (OC, c, and M? 0 
and further all values are given 


long the streamlines 
are randonm constants 
dimensionless form: 


is written 
) . Here 
enthalpy h 


♦Numbers in the margin indicate pagination in the foreign text. 
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relates to heat content h^, values with velocity dimensionality 
relate to the maximum velocity of* steady outflow in a vacuum 
and pressure relates to total pressure p 0 . Under these 


conditions, stream function ^ should be related to 
where 1 is the characteristic linear dimension. 



Flow (1.1) is converted to plane (u, v) of a travel-time 
curve and physical plane (x, y) with equations: 



( 1 . 2 ) 



«/- 

' x p 



(1.3) 



In general, it is impossible to make flow (1.1) adjacent /73 
to the even flow by means of a curved shock because the solution 
is not arbitrary. We will discuss the possibility of joining 
(1.1) with an even flow through an intermediate region of uneven 
flow (Figure 1). If this is set up as an exact problem, it 
can only be solved numerically. To construct an analytic, 
approximate solution we assume that leading shock AC (See Figure 
1) is weak and the flow in fii is isoentropic. In , we have a 
vortex flow. This flow is described with equation (1.1) and is 
joined with fix by means of curved shock wave AB , which diffuses 
in the uneven flow. 



Figure 1 
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In variables h and M>, wave AB (see figure 1) can be 
represented as h = H(\l?) * From (1.3), we have the following in 
the physical plane: 


.v = 


(>l 5 n — t) 


I— a -?-■ a- : .< v:- c F‘ V 

— ,H*r\ /(//) = .v =‘ Fdh,% 

• . *••• .; l •• • -r % • 


(1.4) 




2 (x — I) hlv-t) 


l—o 


’■ h:*-' = urn.; 


0-‘ 


We use ^ to designate the square of the modulus of velocity in 
5 0, Si, 2 are the slope of shock wave AB and velocity vectors 
fli and 9.2 to the symmetry axis. Allowing for the fact that the 
flow in fii is isoentropic, we can then write the laws of momentum 
and energy conservation on wave (1.4) as: 


; l c°s (P — ft,) 1 *; 

. tg (B — fl ) — x ~ 1 HJSLZ-ID. -sin-tP — A.) , I. 

FF F ^ ' x + ! sin (P—0,) cos (P - O a ) ; > * • 

i (1 - l)F = c (*. - [l. - tg (P-^0>) sin (O; -— Ox)]* 


(1.5) 


The law of conservation of mass is satisfied, since the 
coordinate of ^ — stream function — is continuous during passage 
through the shock. The ratio tg/9 = y’^/x\ gives the equation 
for the wave 



( 1 . 6 ) 




The set of equations, (1.1), (1.5), and (1.6), and ratio 




VT=rii ’ 


(1.7) 


obtained from (1.2) are a closed system of equations for z, H, $ , 
Si, S 2 and %. Solving this system defines flow in 9z , the 
position of shock (1.4) and the distribution of gas dynamics 
parameters on the wave front. 
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2. To construct an approximate solution, we must make /74 

the transition to travel-time curve variables in the gas flow 
equations. For problems with axial symmetry, this transition is 
usually useless when studying exact problems and methods, since 
the equations become nonlinear. However, this transition can be 
quite convenient when constructing and analyzing approximate 
solutions . 

In the statement for flow isoentropicity , i.e. h = h(p), 
taken from the Eiler equation transformed to variables % - u 2 + 

+ v 2 = 1 - h, n = v 2 = (1 - h)sin 2 -3i, it follows that there is a 
function 0(£,T|) such that: 

^W v “2/FTV Jyfe- (2.1) 

Any 0 selected will satisfy the dynamic equations^ identically, 
and from the equation for continuity we obtain the ratios which 
determine coordinate x: 


... _ inr 2 (£ ~ ^ \vn ±J&t ± - %v „•, ,/r 2(£- tp ?♦,, 

. . P . Y 8?1 (« - n ) + v n ) ’ ' Y <p n )’ 

_ • . .V: . * r &%■*>:• ' L- v ;.7’; :r> . 


By excluding X from (2.2) with cross-differentiation, we obtain 
the equation which contains only one function for 0: 

. 2 (£ tj) n(q>*Ti — “ l(2£ — q) -f 2^<p^] <p^ — 2r) (2cp^ + 2cp n ) <p^ — 

— n(2 1 1 — j(l — H)] (<Ps -r Tn)} Tnn ~ 2 1/ (<Pj +. Tn)* .+( 1 ~ n/)<Pn(<PS + T<i) «*0. ^ . 


This is a Monzh-Ampere type equation with a quasi-linear part. 

All terms in the equation contain derivatives from the unknown 

function as products of two partial derivatives, and coefficients 

for the products are polynomials according to variable T| to no 

more than the second power. These properties of equation (2.3) 

make it possible to search for the solution as a series according' 

to the powers of T) . We enter the parameter £ = max f? 

ft, 
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and let 


> 9 ==; 2 ® k 9k (I. n>V, n = eri;- 

' ' 0 


(2.4) 


Inserting (2.4) in (2.3) gives the following equation for 0^: 


X (2|/ 1) 9o,tj ^9/i.Tiri — — 9/i. nj — F a— i- 


(2.5) 


Function F n is expressed through 0^ and their derivatives are 
expressed with k <_ n . 

F n — (2s/ — 1) (9«.n ” T l < P*.nn) 9n— »,n ~ 2^Tj9 < ,{,r l 9n— «,«n "r 2 t) 1(1 7 — 2 S/).9«.fc 4“ 
i+n/^#.n+ 2 6tl9«.ul9n-*.iinT- l 2 S9«.U~ 2 t l?*.W — n/ 9 «.ti— (1“ 4 S/) 9*.l) 9»-».n+- 
S+ 2n9«,tiv9n-#.nl + (219^1 + 2l f 9*. W9n-».« 

*■ 

Here, repetitive index s indicates summation from 0 to n. If 
2%j - 1^0 (i.e. Mach number 1), formula (2.5) is a 

recurrent ratio for 0^. Obviously F_^ = 0. Then, from (2.5), 
0o = A 2 x r| 2 /2 + B, where A and B are random functions of 
According to (2.1), the flow corresponding to this solution is 
written with the following formulas: 


x — 


^-(2lj~\)dl if = 2M"-iy 




These formulas determine flow in a channel with variable cross- /75 
section area in a one-dimensional approximation (isobars for % - 
= const are the straight lines for x = const). Further 

. 9 * == -4'*) (5 2*/) ^4^' 4- 3M*J -3l f ' 

The flow corresponding to approximate solution 0 = 0o + £0i is 
described with the following equations: '> 
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( 2 / 6 ) 


7 x tl -f- J (2s/ — ljds.;,: 

' * :•* &• . • j,/.,' ?♦* ; .V. 

;/• 1=5 2 M 8 n t ll^s { aa * t a '*) + 4 ( i ■+ g/j 44' .+ 3/4*] tj* + 

•-V • r ..+‘2(5' + 


Hfere £ is set equal to one, which actually means a return to 
undeformed variable T] . According to (2.6), the isobars are 
second order curves. 


For an rough estimate of approximation process accuracy, we 
compare it with certain exact solutions for gas flow equations. 
Since the basic approximation is an expansion according to the 
angle of flow divergence, a solution in which angle remains 
unchanged in the trajectories (i.e. the flow frotfkthe gas source. 
This flow, by the way, has the simplest structure) is the natural 
choice for comparison. 



Figure 2 

Figure 2 shows graphs of the approximate solutions’ relative 
deviation <S(&) from the exact solution as a function of the angle 
of cone G for two Mach numbers at the inlet. Solid lines (x = 
const) are a one-dimensional approximation and dotted lines are 
approximation (2.6). As Figure 2 shows, the near-axial region is 
described most exactly, the second approximation nearly doubles 
accuracy, and when the Mach number at the inlet (and consequently 


/Jr 
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throughout the entire flow) increases, error decreases. 

3. To construct the flow in fii we use approximate solution 
(2*6). Wave AC (Figure 1) is assumed to be weak, and the angles 
of turn on it are small. It is precisely in this case that 
expansion (2.4) is valid. In variables % and t] , weak shock AC is 
described with the equation: 



From the laws of conservation of momentum and mass and the laws 
for joining the flow on AB, we obtain: 




* 


4*(n.ui 



1AC) J 


(3.2) 


B" _ 2 H ~ 1 + J id b 2, C> p + T l/ c . 2 ryf 

# L 2 V ’Ixc - *Uc A \ 21/ 1 v- ~ 2x - X) 

= ( x T "■ nic - nic, v = 2^.- ^ r^ c ) f 


where y AB and T 1 AB are functions of determinable with /76 

system (1.1) and (1.5) - (1.7); index <» relates to the 
corresponding leading undisturbed flow parameters. 

The general system of equations (1.1), (1.5) - (1.7), (3.1) 

and (3.2) determines total flow, i.e. gives the solution to the 
problem of joining vortex flow (1.1) with even flow. 

Shock wave AB reaches the axis if y = 0 in (1.3), i.e. if 
either M? = M?o (a) or z* - 0 (b). In first case (a), we have po = 
= 0 for total pressure on the axis, which is impossible. We will 


9 


0 . 


consider the second variation, z* = 0. For (1.1), this is a 
singular point through which the single-parameter family of 
equations z = “ ho + k(ho — h) 2 passes (ho is the value for 

enthalpy on the axis and k is the parameter). Here, according to 
(1.2), v = 0 on the axis. It follows that the shock wave reaches 
the axis and approaches it at a straight angle. 


Calculations are made from the axis. Given M m and (X, we find 
parameter k, which fixes the integral curve in the family of 
solutions passing through the antinode point. Further, departing 
from the singular point according to the analytic solution, we 
solve the general system of common differential and 
transcendental equations (1.1), (1.5) - (1.7), and (3.1) - (3.2). 

At each integration step for the differential equations, the 
system of transcendental equations (1.5) is solve** by iteration 
[2]. Here a solution from the preceding step is selected as the 
initial approximation for the iteration process. Calculation 
continues up to y = y max , whose value corresponds to flow lines 
\ax = 


4. The intensity of shock wave AB decreases as y^g increases 
and falls to zero when y = yj„ ax . Intensity changes due to solid 
wall distortion and uneven flow in ih . At values of <X <_ 1, the 
unevenness of flow in is insignificant (when CX = 0.1, the 
maximum value of Atfi/A^, on the external side of the shock is 
0.06, and when (X = 0.01, A>J i / A^j, = 10~^) and shock distortion is 
due basically to wall distortion. In this case leading wave AC 
is almost a characteristic curve and the zone has significant 
dimensions. As (X increases, the flow in fii becomes more uneven, 
AC begins to differ from the characteristic curve, and zone fti 
decreases. Here the vorticity in Qz can be explained not only by 
wall distortion, but also by uneven flow in front of the shock. 
Note that here the hypotheses on which the solution in are 
based are met with decreasing accuracy. 


10 


i 

I 



Figure 3 


We use the solution obtained to describe the flow in 
airscoops constructed by an excision along the streamlines of an 
axisymmetric flow [3, 4]. We will discuss a sector excision with 
angle 20 (Figure 3 a). The combined loss coefficient for total 
pressure is determined according to the formula 

|«i= SJpoJS, f. (4.1) 
where S is the area of the capturable gas jet. 

The average recovery factor for total pressure for this 
diffusor formula is 



(4.2) 


Figure 3 b shows ]1 as a function of the angle of sector 0 at /77 
different values for S, and Figure 3 c shows JJ as a function of 
index C(. The influence of 0( on the recovery factor increases 
with 0. In Figure 4 the values of 0^ are determined with the 
equation 0j = S^/y max , which shows that y <_ y max . 


* 
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To evaluate the level of compression for the gas jet entering 
the airscoop, we determine pressure on AD (see Figure 1). On AD 
h - const = h j ♦ Then 


We obtain the following for total compression 

... . W . n ‘- ,'y 

Ci'== 

where JJ X is determined in (4.1). 



The characteristics of intake diffusors with other excision 
shapes can be determined in a similar manner. Figure 4 shows a 
diffusor with a rectangular excision which can be characterized 
with two parameters: the area of captured jet S and angle 0, 
shown on Figure 4 a. Figure 4 b shows the curve for the ratio of 
diffusor recovery factor (Figure 4 a) to the analogous value / 78 
(4.2) characterizing the influence of excision shape on 
aerodynamic characteristics. 




12 


REFERENCES 


Dulov, V. G., Shchepanovskiy, V. A., Ob odnom tochnom 
reshenii uravneniy osesirametrichnykh neizoentropicheskikh 
techeniy [On On© Exact Solution for Equations for 
Axisymmetric Anisoentropic Flows], Izv. Sib. otd. AN SSSR, 
Seriya tekhnicheskikh nauk, No. 13, Vol. 3, 1976. 

Berezin, S. I., Zhidkov, N. P., Metodv vyehisleniva [Methods 
of Calculation], Vol. 2, Moscow, "Nauka" Press, 1966. 

Shchepanovskiy, V. A., Postroyeniye pros t rans t vennogo potoka 
iz ploskikh techeniy [Constructing a Three-Dimensional Flow 
from One Dimensional Flows], Isv. Sib. otd. AN SSSR, Seriya 
tekhn. nauk, No. 8, Vol. 2, 1972. 

Gutov, B. I., Zatoloka, V. V ., Convergent Intake Diffusers 
with an Initial Shock Wave and Additional External 
Compression, in: Aerof izicheskive issledovaniva 

[Aerophysical Research], Vol. 2, Novos ib irsk 1973 . 



